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SUMMARY 

The present paper concentrates on the derivation emd intepretatlon of the 
necessary conditions that a sailplane cross-country flight has to satisfy to 
achieve the maximum global flight speed. Simple rules ace obtained for two 
specific meteorological models. The first one uses concentrated lifts of vari- 
ous strengths cmd unequal distance. The second one takes into account finite, 
non-uniform space amplitudes for the lifts and allows, therefore, for dolphin- 
style flight. In both models, altitude constraints consisting of upper and 
lower limits are shown to be essential to model realistic problems. Ntanerical 
examples illustrate the difference with existing techniques based on local 
optimality conditions. 


INTRODUCTION 

The problems associated with the optimization of sailplane flight paths to 
achieve maximum cross-country speeds have recently received special attention 
in the literature. This has been stimulated by the modern competitive soaring 
which consists almost exclusively in racing and by the development of high 
performance sailplanes allowing fcx* new, highly efficient flight techniques. 
Starting with the now classical MacCready [l] results, most of the investiga- 
tions have be^ concerned essentially with local optimization problems, that 
is, finding the optimum flight strategy for various specific situations 
encountered in a short section of a flight [l to 10]. 

In recent papers (2, 4, 5, 8) the optimum speeds to fly in a variety 
of atmospheric vertical velocity distributions have be«i determined from the 
basic assumption that the corresponding flight segments had to be crossed 
with zero net altitude loss. Conditions under which a transition from the 
circling mode of climb to the dolphin or easing modes has to be decided have 
been examined [ 4 ]. Although such results yield extremely valuable guidelines 
for selecting a flight strategy, they only optimize the speed over a limited 
portion of the total flight. 

It is well known, however, in optimization theory that a succession of 
locally optimisn solutions does not, in general, lead to a globally optimum 
result [ill. It is worth pointing out that a globally optimum flight strategy 
can only be determined if the assumption is made that the distribution of 
atmospheric velocities over the whole flight path is known in advance and 
is independent of time. Although this is never achieved in practice, it is 
felt that the derivation of global optimality conditions allows for a new 
insight into the sailplane flight technique by giving a posteriori the deci- 
sions that the pilot should have taken and the influence of factors that have 
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been up to now neglected in the analysis, such as the effect of the unequal 
distribution and strength of the lifts euid the effect of minimum and maximum 
altitude limitations. Such altitude constraints reveal an essential ingredient 
in the formulation. Their necessity appears as follows. If they are absent 
and if the lifts are of unequal strength, the globally optimum solution turns 
out to be trivial and consists of a glide intil the maximum lift existing 
along the pa^h is reached and a climb to an altitude that allows completion 
of the task, the speed on both segments corresponding to the MacCready setting 
for that strongest lift [l2l. 

The present paper provides simple rules for global optimality for two 
simple atmospheric models. These appear to be in agreement with the techniques 
intuitively used by good competition pilots. 


PROBLEM FORMULATION 

In both atmospheric models used in the following, the horizontal (wind) 
velocity of the ait mass is assumed to be either zero or to be taken into 
account by an appropriate modification of the polar equation. The vertical 
velocity (lift) of the air mass cj is defined by the so-called netto value. 
It is constant in the vertical direction between the altitude limits. The 
flight path is supposed to be constituted by a succession of segments of vari- 
able lengths in which the air mass exhibits vertical velocities c^ which are 
constant along a given segment but vary from one to the other. The altitude 
constraints consist in constant upper and lower limits denoted h and h. 

Note that variable altitude limits could be easily Incorporated. For simplic- 
ity the lower altitude limit is taken as zero (h > 0) . 


Concentrated Lift Model 

In a first model, the lengths of the segments where a positive vertical 
velocity is encountered are supposed to be negligible, that is, the lifts are 
considered as concentrated. The air mass between the lifts is supposed to be 
stable. Climbing is therefore achieved only by circling at fixed locations 
corresponding to the lifts. If a lift is not used, its crossing Is supposed not 
to affect the glide and the dynamical aspects of the transition from gliding 
to circling are not considered. The model is illustrated in figure 1 which is 
drawn in the vertical plane. The signs of the velocities are taken according 
to the positive sign of the axes. 

The flight consists of a succession of climbs in the selected lifts fol- 
lowed by a glide at constant speed which possibly crosses discarded lifts. The 
pilot controls the selection of the lifts where he decides to gain altitude, 
the amount of altitude gained, and the speed to fly between the selected lifts. 
During the climbs the vertical speed of the sailplane is simply taken as the 
algebraic sum of the air mass vertical velocity and the minimum sinking speed 
of the glider. The increase of sinking speed due to variations in bank angle 
is not considered explicitely and should be incorporated in the polar definition. 
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The variations in altitudes are given by 


^21+1 - ^2i = Ahi 

W(Vi) 

^2i+2 - ’^2i+l = ^-i 

Vi 


i = 0,l,...,n-1 (1) 


where the sinking speed w(V;^) is given by the polar equation. The classical 
quadratic approximation has teen used for numerical examples 

w = A v^ + B V + C (2) 


The time spent at each step consists of the sum of the time used in climb- 
ing and the transition time between lifts 


Ahi 

t2i = — 

ci + «m 

^i 

t2i+1 = — 

Vi 


i = 0,1,...,n-1 (3) 


The achieved rate of climb is the sum of the air mass velocity C£ and the 
glider minimum sinking speed The constraints on the altitude and altitude 

gains are expressed by 


A control constraint 


Ahi = 0 


(4) 


Initial and terminal 
constraints, say 


ho = 0 h2n = 0 


(5) 


Altitude constraints 
at each step 


h2i+i ^ h; h2i - 0 


( 6 ) 


1 ~ 1 


The mathematical problem can be treated by the classical discrete optimal con- 
trol theory (ill and consists of finding the sequence (s) Aho, Vq, Ah^ , v^ , 

Vn_i which satisfy the relations (1), (4), (5), and (6) and minimizes the total 
flight time 


T 


n-1 

y 

i=0 


'^2i * 


(7) 
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Distributed Lift Model 


In this second model, the length of the segments is always finite and non- 
vanishing. Between lifts the air m 2 uss may present negative vertical velocities 
so that <uiy desired air mass vertical balance can be achieved. The model is 
illustrated in figure 2. An important difference with the preceding model is 
brought by the possibility of crossing a lifting segment at a horizontal speed 
Vj^ less than the speed corresponding to the minimum sinking speed w,p. This 
is achieved by using the equivalent polar illustrated in figure 3 and already 
used in other similar works [ 4 I . For horizontal speeds less than that corre- 
sponding to the minimum sinking speed, the sinking speed remains constant. This 
appears to be sufficiently accurate to simulate the transition from pure dolphin 
flight to essing or circling or a combination of equivalent manoeuvers achieved 
to cross a lifting area in the time regjired to gain a certain amount of alti- 
tude. Note that the same approximation is used as the basis for speed control 
in some modern instrumentation. In the numerical applications, the quadratic 
approximation (2) remains applicable at speeds higher than v(w„|) . 

The variations in altitude are governed by 

W(Vi) + Ci 

hj^+l — hj — ' ' i i — 0,1,...,n— 1 (8) 


If Vi ^ v(w,n) , then w = Wni = Constant and thus the altitude gain is entirely 
controlled ty the equivalent horizontal speed V£ . The altitude gain no longer 
appears eis an explicit control variable. The time spent in each segment is 
given by 

^i ^i+T - hi 

ti = — = (9) 

Vi W(Vi) + Ci 

while the altitude constraints read 


At initial and terminal hQ = 0 h„ = 0 (10) 

points 

At each step 0 ^ hi S h i = 1,2,...,n-l (11) 


The mathematical problem consists of finding the sequence(s) Vq, vt , ...» 
Vn_T satisfying the relations (8), (10), and (11) and minimizing the total 
flight time 


T 


n-1 



i-0 


( 12 ) 
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NECESSARY OPTIMALITY CONDITIONS 


The first-order necessary conditions for optimality can be deduced by the 
classical methods of discrete optimal control [ill . Such methods have been 
used in previous work [4, 6] . A detailed treatment can be found in [l2, 13] 
for the two atmospheric models presented here, as well as for certain more com- 
plex situations. The conclusions are summarized as follows. 


Concentrated Lift Model 
The Hamiltonian turns out to be 

&hi 

H(h2i, Ahi, P21+1, P°/ 2i) = p° + P2i+i Ah^ 

ci + 

*-i W(Vi) 

H(h2i+i, Vi, P2i+2» 2i+l) = pO _ + p2i+2 i'i 

Vi Vi 


(13) 


It has to be maximized with respect to Ahi or vi for each i. The so-called 
adjoint variables Pi have to satisfy the relations 


pO - Pl =0 

P2i+1 ■ P2i+2 = ^2i+1 
^2i+1 (1'2i+l - h) = 0 
P2i ■ P2i+1 = "^2i 
^2i ^2i = 0 


pO S 0 

i = 0 , 1 , . . . , n- 1 
X2i+2 = 0 (14) 

i = 1 ,2, . . .,n-l 
X2i ^ 0 


where the Xi are Lagrange multipliers. Prom those conditions, it can be 
shown [l2] that a reduced set of adjoint variables t)q, Dif , hn-v 'I^o' '^’1' 
..., 4'n-l Lagrange multipliers Pi, P2» ...» Pn-i can be derived which, 

in an optimal solution, have to satisfy the following conditions 


H2i+1 (^>2i+1 - h) = 0 P2i+1 - ® 

^2i ^'2i = 0 V>2i = 0 

Ahi = ® hi ^ 0 hi Ahi ~ 

hi+1 = Hi + Ci - ci+1 + P2i+i - M2i+2 
dw(vi) 

Vi - w(vi) + il^i = 0 t^i > 0 

dVi 


0,1,..., n-1 

(15) 

1,2,...,n-1 

(16) 

0,1 ,..., n— 1 

(17) 

0,1,..., n-2 

(18) 

0,1,..., n-1 

(19) 
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't’i+1 = 'i'i - W2i+2 + ^2i+3 

i = 0,1,... ,n-2 

(20) 

il^i = Hi + 9i + w„ + M2i+1 

i = 0 , 1 , . . . , n- 1 

(21) 

Distributed Lift 

Model 


The optimal solution (hi, Vi) must be such that the Hamiltonian 


^i 

H(hi,Vi, Pi+T, pO, i) = pO _ + 

n 

w(Vi) + Ci 

« . 

(22) 

^i 

Vi 


i = 0,l,...,n-l 


maximized with respect to vi for each i. 
satisfy the relations 

The adjoint variables pi 

have 

0 

II 

r— 

a 

0 

a 

0 

VII 

0 

a 


.H) .(2) 

Pi “ Pi+7 = -Xi + Xi 

i = 1,2,...,n-1 

(23) 

Xi hi = 0 

0 

VII 


, f2) - „ 

Xi (hi - h) = 0 

2 

Xi ^ 0 



where and Xj^ are Lagrange multipliers. Again a set of reduced vari- 

ables and Lagrange multipliers 

/. , 0) H) (2) (2)\ 

('pQ» • • • » • • • fPn-l 1 


allows to present the optimality conditions in a more suitable form which turns 
out to be [ 1 3 ] 


( 1 ) 

Wi hi = 0 


( 2 ) 

Ui (hi - h) = 0 


( 1 ) 

‘i 

( 2 ) ^ 

Mi ^ 0 


i = 1 , 2, . . . ,n-l (24) 


(Vi) 

Vi w(vi) 

dVi 


- Ci + 'I’i » 0 


i = 0, J , . . . ,n-1 (25) 


360 


(26) 


'I'i = - Mi + Mi 1 = l,2,...,n-l 

'l>i > 0 'J'i “ Wfli + Ci i = (27) 


PHYSICAL INTERPRETATION 


The interpretation of the two sets of optimality conditions follows similar 
lines. A first conclusion is drawn from equation (19) or (25) which governs 
the speed to fly in a segment. Note that equation (25) reduces to equation (19) 
if the air mass vertical velocity Ci is zero, as assumed in the first model. 
From figure 3, the reduced adjoint variable appears to correspond to a 

classical MacCready setting and indeed equation (25) appears in most other works 
on optimization [4, 5, 8]. In the following, the notation MCS(Ci) denotes the 
setting corresponding to an air mass velocity Ci as defined by equation (25) . 
The next interpretation concerns the Lagrange multipliers Mj in equation (15), 


(16), or (24). 

ing ( V2i+1 or 
altitude limit) 


These multipliers are zero in entering 
( 2 ) 


(^2i 


( 1 ) 


or Mi 


or leav- 


a segment if the LAL (lower altitude limit) or UAL (upper 


is not reached. 


From equation (20) or (26) the important conclusion is drawn that the MCS 
cannot change from a segment to the next one unless either the UAL or the LAL 
has been reached, that is, if one of the Lagrange multipliers m becomes nega- 
tive. If, and only if, the UAL is reached, then the MCS may be reduced. Con- 
versely, the MCS may be increased only if the LAL is touched. To proceed far- 
ther with the interpretation requires distinguishing between the two models. 


Concentrated Lift Model 

The reduced adjoint variables appears from equation (17) as indicators 

of whether the lifts c^ may be used (rii = 0) to climb or not (hi > 0). With 
these results in mind, it becomes easy to deduce from equations (18), (20), 
and (21) the logic for deriving iteratively the optimum solution. 

Consider the beginning of the flight at a location where a lift A exists. 
Denote by B the first lift stronger than A along the flight path. The follow- 
ing iterative reasoning can be made. If B can be reached from the UAL in A with 
a MCS(A), that is, a MCS corresponding to the lift A, then one has to climb 

in A just enough to reach B at LAL with a MCS(A). If B cannot be reached with 

MCS (A) even when climbing at UAL in A, then climb to UAL in A and look fo' the 
next best lift, denoted C, between A and B. Evidently C < B and C S A. Try 
to reach B at LAL with a MCS inferior to MCS(A) but superior to MCS(C). If this 
is impossible because C cannot be reached, then restart the reasoning with A 
unchanged and B replaced by C; if this is impossible because B cannot be reached, 

then take a MCS(C) to reach C and restart the reasoning with B unchanged and 

A replaced by C. 
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Consider now the case where a lift A is stronger than all remaining ones 
on the flight path. Denote by B the strongest of the remaining lifts. Unless 
in the special case where the task could be ended from A with a MCS(A) without 
climbing up to the UAL, one necessarily has to climb up to the UAL and take a 
MCS superior to the MCS(B) but equal or inferior to the MCS(A). If the task 
cannot be ended, reach B with a MCS(B) and climb in it up to the UAL. Repeat 
the reasoning in B for the next strongest lift. If B cannot be reached from 
A with MCS(B), look for the best lift between A and B, denoted C, and try if 
B can be reached at LAL with a MCS between MCS(B) and MCS(C). If necessary 
climb in C if the MCS is equal MCS(C). If it is impossible to reach C at LAL 
with a MCS(C) look for the best lift, say D, between A and C and repeat the 
reasoning. If there is no lift, tryMCS(O). If it still does not work, the 
flight is evidently impossible. 

A combination of the two reasonings proposed above for increasing or 
decreasing lifts allows the construction of the optimal solution. Note that 
it is not necessarily unique. 

It is worth pointing out that the optimal solution leads always when going 
from a lift A to a lift B to use a MCS corresponding to the weaker of the two 
lifts. If A > B the UAL has to be taken in A while B has to be reached at LAL 
if A < B. Similar conclusions have been obtained independently by a variational 
approach in [ 1 4] . 


Distributed Lift Model 

In this model, the decision to gain altitude in a lift is dictated by equa- 
tion (27). If ^ w^ + Cj^ the speed vj is larger than (w,p) and there- 
fore is uniquely determined by equation (25) . The lift has then to be crossed 
at the corresponding speed. This appears to be a pure dolphin mode. If, and 
only if, 'i’i = w„, + cj the decision of climbing may be taken as the speed 
becomes equal to or smaller than the speed v(Wjp) . Indeed, due to the form used 
for the equivalent polar, the speed vj cannot be computed by equation (25). 
Its value is dictated by the need to gain a certain amount of altitude in that 
lift, given by 


1 

Ahi = ti = (w„, + Ci) — (28) 

Vi 

If the speed is inferior to vj (Wj„) , Ahj^ is larger than the value 

obtained in pure dolphin mode which implies that some sort of manoeuver like 
essing and/or circling is achieved while flying through the segment. 

Except for the process of gaining altitude in a non-dolphin mode, the con- 
clusions reached for the preceding model are still valid. From equation (26) 
the MCS may not be changed unless the altitude limits are reached. It may 
increase only if the LAL is touched and may decrease only when reaching the UAL. 
The process for constructing numerically an optimum solution is as follows. 
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start by trying to use the MCS of the best existing lift, say A, for all 
the segments. If the LAL is not reached, increase the setting until either 
the task can be ended or the LAL is reached. At that point, the MCS may be 
increased. 

If the MCS corresponding to the best lift allows reaching the LAL before 
the segment where it occurs, look for the strongest lift between the present 
point and A. Denote it by B. Then try to reach A at the LAL with a MCS(B) . 

If this is not possible, climbing in B is allowed. If B cannot be reached with 
a MCS(B), look for the best lift between the present point and B and repeat the 
reasoning as necessary, keeping in mind the rules that allow climbing in a lift 
and those that allow changing the MCS. 


NUMERICAL EXAMPLES 
Concentrated Lift Model 

As a simple example, a 300 km flight is schematized in figure 4. The lifts 
are equidistant (10 km) for simplicity although it is by no means implied in 
the preceding rules for optimality. The lift strengths are indicated in m/sec 
along the y-axis. They increase progressively during the flight, then decrease, 
but cire in general unequal. The altitude limits are 0 and 1000 m. The sail- 
plane polar is approximated by 


w = -1.65 10-3 v2 + 61.6 10-3 ^ - 1.026 (m/sec) 


and corresponds to a dry open class ship. The optimal strategy for that lift 
distribution and altitude constraints is illustrated in figure 4. The MCS for 
each glide is indicated. It follows as a simple and systematic application of 
the rules for optimality established above. Note that the flight strategy con- 
sists in hitting systematically the altitude limits, except at 110 km and 170 km 
where the altitude just necessary for reaching the next best lift at LAL is 
gained. Note also that the MCS is not always equal to the strength of one of 
the two lifts defining the glide. Finally, note that this example justifies 
the practical rule of flying "low" when the lifts are improving and keeping 
"high" when they are deteriorating. The cruising speed for that flight is 
81.01 km/h. 

As a test for the sensitivity of the speed with respect to the MCS, the 
same problem has been solved with the additional constraints of keeping the same 
MCS which implies the same speed between any two selected lifts. The optimal 
MCS for that situation has been obtained in fl3l in the form 

L 

i]) = (29) 

Y s-i 

^ Ci + Wm 
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where L is the total length of the flight (300 km) and cj + w„i is the 
achieved rate of climb in each selected lift. Note that the selection of the 
lifts is highly dependent upon the altitude limits. In the present example 
>1^ = 1.23 m/sec and the corresponding constant speed between the thermals is 
133 km/h. The new constrained optimal speed becomes 79.51 km/h which differs 
by only 1.85% from the exact optimum. Although some restrictions have to be 
mentioned (see [l3]) concerning the applicability of equation (29) in a general 
case, it indicates clearly that in an atmosphere corresponding to the present 
model, the MCS is much less important than the correct selection of the lifts. 
This is again well known for many competition pilots [ 1 5l . 


Distributed Lift Model 
The flight polar has been approximated by 


w(v) = -1.896 10"3 v2 + 77.8 10“3 v - 1.27 (m/sec) 


which corresponds also to a dry open class glider and to the model used in [4l 
and [5] . Three distributions of lifts have been selected iind are presented in 
figures 5, 6, and 7 and denoted flights I, II, and III. These flights are all 
200 km long and correspond to different weather conditions. In flight I the 
lifts are relatively concentrated except at two places and their strengths are 
rather different from each other. The length of the lifting zones represents 
36% of the total which is rather critical for the transition from thermaling 
to dolphining [4]. The air mass balance is positive, that is, the average over 
the distance of the air mass (netto) vertical velocities yields 0.39 m/sec. 

In flight II the lifting zones represent 31% of the distance and their strengths 
are much more similar to each other. The lifting and sinking areas exactly bal- 
ance each other? that is, not only the average vertical velocity is zero, but 
the air mass is organized in a succession of cells which are 20 to 40 km long 

where the exact air mass balance is also achieved. This allows for using the 

classical rules for local optimality ( 4l in crossing these cells and compares 
with the globally optimal solution. In flight III the lifts are weaker and 
their strengths still closer to each other. The lifting zone represents 49% of 
the total. The air mass balance yields 0.236 m/sec and the lifting zones are 

again organized in cells in which approximately the same air mass balance is 

maintained. For each of these atmospheric models three upper altitude limits 
have been considered h = 1000 m, 1500 m, and 2000 m. The LAL has been kept 
at h = 0 which is evidently not necessarily the ground level. 

The numerically obtained optimal solutions ate illustrated in tables I, 

II, and III in digital form and in figures 8, 9, and 10 in graphical form. The 
satisfaction of the optimality conditions described above are easily verified. 
The lifts in which gaining altitude in circling or essing are indicated as well 
as the corresponding equivalent horizontal speed which is then smaller than 
the speed of minimum sink vCw^j) = 20.52 m/sec. In the other segments, crossed 
in dolphin mode, the optimum MCS is given. Note that = 0.47 m/sec. The 
influence of the altitude limits is illustrated by the following table: 
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Ah = 1000 m 

Ah = 1 500 m 

1 

Ah = 2000 m 

Ah = Unlimited 

Flight I 

V = 94.5 

97.94 

100.19 

100.57 

Flight II 

V = 73.76 

81 .2 

83.10 

84.20 

Flight III 

V = 85.87 

87.98 

88.16 

88.16 


where Ah is the allowed altitude range and v is the optimum average speed 
in km/h. The application of various non-globally optimal flight strategies 
based on the use of existing rules for optimizing the speed in each individual 
meteorological cell [ 4 , 5l resulted in average speed frau 5 to 15% inferior 
depending on the allowed altitude limits and on the various conditions selected 
in applying these rules . 


OONCLUSIONS 

Simple rules for obtaining numerically the optimum flight strategy in two 
meteorological models have been obtained. Their applications reveal that the 
altitude limits imposed for the flight may have, as known from experience, a 
much more significant influence on the achieved speed than the selection of MCS. 
Additional investigation is required to determine the relation beween the vari- 
ous possible weather profiles and the optimum flight strategies as well as 
altitude limits influence. The in flight recording of such atmospheric profiles 
over rather long distances would allow for studying systematically the optimum 
solution corresponding to a number of classical situations. 
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SYMBOLS 


A,B,C 


hi 

h 

H 

Pi 

ti 

T 

Vi 

w 


^ i »11 i»U i 

<^i 

i 

UAL, LAL 
MCS 


flight polar coefficients 

air mass netto vertical velocity 

altitude 

upper altitude limit 

Hamiltonian 

flight segment length 

adjoint variable 

elapsed time 

total flight time 

horizontal speed 

sailplane sinking speed 

minimum sinking speed 

Lagrange multipliers 

reduced adjoint variable = MCS 

index of a flight segment 

upper (lower) altitude limits 

MacCready setting 
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TABLE I.- OPTIMAL SOLUTION FOR FLIGHT I 



H 

- 1000 m 


H 

- 2000 m 



max 



1 max 



Segment no . 

MCS 

out 

B 

Mode 

MCS 

out 

B 


1 

0.53 

425 

0.6 

C 

0.53 

425 


C 

2 

0.53 

0 

30. 8 

D 

0.53 

0 


D 

3 

1.03 

481 

1.0 

C 

1.03 

481 


C 

4 

1.03 

0 

34.8 

D 

1.03 

0 


D 

5 

2.03 

732 

2.7 

C 

2.03 

732 

2.7 

C 

6 

2.03 

288 

41.7 

D 

2.03 

288 

41.7 

D 

7 

2.03 

213 

38.4 

D 

2.03 

213 

38.4 

D 

8 

2,03 

233 

34.8 

D 

2.03 

233 

34.8 

D 

9 

2.03 

0 

44.7 

D 

2.03 

0 

44.7 

D 

10 

3.03 

164 

34.8 

D 

3.03 

164 

34 . 8 

D 

11 

3.03 

1000 

18.1 

C 

3.03 

1781 

9.3 

C 

12 

1.57 

646 

47.8 

D 

3.03 

1396 

55.3 

D 

13 

1.57 

356 

45.0 

D 

3.03 

1069 

52.8 

D 

14 

1.57 

249 

35.0 

D 

3.03 

826 

44.7 

D 

15 

1.57 

350 

31.0 

D 

3.03 

749 

41.7 

D 

16 

1.57 

712 

21.0 

D 

3.03 

913 

34.8 

D 

17 

1.57 

53 

42.0 

D 

3.03 

121 

50. 3 

D 

18 

1.57 

0 

35.1 

D 

3.03 

0 

44.7 

D 

19 

4.61 

268 

42.2 

D 

4.03 

370 

38.4 

D 

20 

4.61 

1000 

27.0 

D 

4.03 

2000 

12.3 

C 

21 

1.38 

428 

43.8 

D 

2.69 

1365 

51 . 1 

D 

22 

1.38 

0 

40.7 

D 

2.69 

858 

48 . 5 

D 

23 

1.53 

858 

8.9 

C 

2.69 

1056 

32.1 

D 

24 

1.53 

437 

38.4 

D 

2.69 

506 

45.7 

D 

25 

1.53 

0 

41.7 

D 

2.69 

0 

48.5 

D 

Speed 


94 . 54 

km/h 


100.19 

km/h 



■ altitude at the end of the segment (meters) 

Mode * D for dolphin 

C for climbing with v. < v . (w ) 

1 I m 

MCS, » m/sec 


368 















FOR FLIGHT II 


H » 2000 m 

max 

MGS 

^out 

B 

Mode 

0.53 

447 

2.9 

C 

0.53 

284 

30.8 

D 

0.53 

541 

20. 5 

D 

0.53 

0 

38.4 

D 

1.03 

617 

8.3 

C 

1.03 

0 

38.4 

D 

2.03 

2000 

5.0 

C 

2.03 

1 398 

47.6 

D 

2 . 03 

1 1 65 

44.7 

D 

2.03 

1314 

38.4 

D 

2.03 

1 149 

30.6 

D 

2.03 

450 

44.7 

D 

2.03 

149 

47.6 

D 

2.03 

2000 

5.5 

C 

1 .53 

17 11 

44.7 

D 

1 .53 

1055 

41.7 

D 

I .53 

1215 

30.8 

D 

1.53 

1 074 

38.4 

D 

1 .53 

200 

41.7 

D 

1.53 

1092 

8.6 

C 

1.53 

0 

41.7 

D 


83.10 

km/h 



the segment (meters) 


< V . (w ) 

T 1 ' m ' 


















TABLE III 


OPTIMAL SOLUTION FOR FLIGHT III 



H 

max 

“ 1 000 m 



H - 2000 m 

max 

Segment no. 

MCS 

ou t 

B 

Mode 

MCS 

^out 

B 

Mode 

1 

0.53 

364 

3. 

6 

C 

0.53 

364 

3.6 

C 

2 

0.53 

200 

30. 

8 

D 

0.53 

200 

30.8 

D 

3 

0.53 

586 

20. 

6 

D 

0.53 

586 

20.6 

D 

4 

0.53 

0 

34. 

8 

D 

0.53 

0 

34.8 

D 

5 

1 .03 

999 

20. 

6 

D 

1.03 

999 

20.6 

D 

6 

1 .03 

594 

47. 

6 

D 

1.03 

594 

47.6 

D 

7 

1 . 03 

37 

41 . 

7 

D 

1.03 

37 

4 1.7 

D 

8 

1 .03 

968 

16. 

6 

C 


968 

16.6 

C 

9 

1 . 03 

4 1 1 

41 . 

7 

D 

1.03 

41 1 

4 1.7 

D 

10 

1.03 

0 

38. 

4 

D 


0 

38.4 

D 

1 I 

1.53 


15. 

3 

C 

1.53 

1 737 

8.8 

C 

1 2 

0.56 

608 

35. 

0 

D 

1.53 

1299 

4 1.7 

D 

13 

0.56 

208 

44. 

9 

D 

1.53 

887 

50.3 

D 

14 

0.56 

587 

20. 

8 

D 


1 046 

30.8 

D 

1 5 

0.56 

0 



D 

1.53 

390 

41.7 

D 

1 6 

1.03 

673 



C 

1.53 

759 

26.2 

D 

1 7 

1 .03 

617 

30. 

8 

D 

1.53 

656 

34.8 

D 

1 8 

1 .03 

0 

38. 

4 



D 

1.53 

0 

4 1.7 

D 

Speed 


8 5.87 

km/h 



88.16 km/h 


= altitude at the end of the segment (meters) 

Mode = D for dolphin 

C for climbing with v. < v . (w ) 

I 1 m 

MCS, = m/sec 


370 























Figure 3.- The equivalent polar. 
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Figure 5.- Lift distribution for flight I. 
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C; (m/sec) 10 1-1 1,5 

2.5 75 10 10 5 




200 km 


11 12 13 14 15 16 17 18 

1.5 -.5 -1 2.5 -1 -.5 1 0 

5 15 5 5 5 15 15 5 


MM M b I 5 I I lb I b I I b I I lifting ZONES = 31 7 o 

Figure 6.- Lift distribution for flight II. 
























DISTRIBUTED LIFT MODEL 
1.5 


liOeBBI 




flight segment no. i 2 3 4 5 6 7 8 9 10 L = ^ li = 200 km 

Cr (m/sec) 10 1 -.5 1.5 -2 -1 1.5 -1 -.5 

li (km) 2.5 7.5 15 15 20 5 10 15 10 10 ^^1 

i — r _ <HH = 9”^?; m 

11 12 13 14 15 16 17 18 ^ sec 


Figure 7.- Lift distribution for flight III. 
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Figure 8.- Optimal i;olutiiin fcr flight I, 






